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The influence of the spin-dependent phase shifts (SDPS) associated to the electronic reflection 
and transmission amplitudes acquired by electrons upon scattering on the potential barrier on the 
Andreev reflection probability of electron and hole excitations for a ferromagnet /isolator/d- wave 
superconductor (FIS) contact and the charge conductance of the FIS contact is studied. Various 
superconductor orientations are considered. It is found that SDPS can suppress the zero-potential 
peak and restore finite-potential peaks in the charge conductance of the F/I/d-wave superconductor 
contact for the (110) orientation of the d-wave superconductor and, on the contrary, can restore 
the zero-potential peak and suppress finite-potential peaks for the {100} orientation of the d-wave 
superconductor. 
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I. INTRODUCTION 

The oscillating character of the spatial dependence 
of the anomalous Green function (GF) in a ferromag- 
net in various hybrid structures containing the ferromag- 
net/superconductor (F/S) interface with a singlet order 
parameter is due to the presence of electron spin sub- 
bands with different values of Fermi momenta in a 
ferromagnetic metal (F) [U, 0, H, 0] (a =Tji) is 
the spin index, which denotes the projection of the elec- 
tron spin on the direction of the magnetic moment of a 
ferromagnet). Such manifestation of the proximity ef- 
fect is the basis for creation of the 7r-Josephson junction 
01 , various spin- valve schemes [1, 0, H, 0, , being the 
main elements of promising superconducting electronics 

The suppression of the Andreev reflection Q| in point 
F/S contacts [1^1 due to the decrease of the number of 
conducting channels is another consequence of the pres- 
ence of spin subbands in a ferromagnetic metal. This 
fact is used to determine the spin polarization of ferro- 
magnetic materials [H, [13, 0, [H, Hyl , to study the order 
parameter s y mmetr y of high-temperature superconduc- 
tors [U, I22I [2p. [241 and to control the spin-polarized 
currents |25l . \2ui . \2T\ . 

Recently an attention was paid to one more property 
of hybrid F/S structures: the influence of spin-dependent 
phase shifts 6^ and 6^ (SDPS) associated to the elec- 
tronic reflection and transmission am plit udes Tq, and da 
on the contact on thermodynamic |28| and transport 
[29l . [sol . [3l| characteristics of hybrid structures with a 
spin-active interface: 



da = ^/D^ expiiOi))] 



'R^ exp{i9a). 



Here D a and i? „ = 1 — D a are transmission and reflec- 
tion coefficients, respectively. Let us note, that parame- 
ters da and are almost insensitive to the appearance 
of a superconducting state and the corresponding changes 
will not be considered in the present paper. 



It has been found that the difference of SDPS due to 
the difference of potential barriers for electrons with dif- 
ferent spin projections a results in the appearance of a 
TT state in the S/FI/S junction (FI is a ferromagnetic iso- 
lator) without taking into account the proximity effect 

The presence of SDPS also leads to the formation of 
spin-dependent Andreev bound states = sign{px,a — 
Px,-a) ^co&{{px,a — Px,-a)/2) in the N/F/s-wave su- 
perconductor contacts (N is a normal metal; Px, a is the 
projection of the Fermi momentum in a ferromagnet on 
the X axis, being perpendicular to the contact plane) 
[Slj . In the tunneling limit these states appear as the 
resonance peaks below the gap in the dependence of the 
ballistic charge conductance on the applied potential V. 

The influence of SDPS on the charge conductance of 
a single-channel quantum point contact of a F/s-wave 
superconductor and that of a multichannel ballistic con- 
tact of a F/I/s- wave superconductor (I is isolator) was 
studied in Refs. [H, [s^l, correspondingly. In Ref. [Ill 
it has been found that for a weakly transparent contact, 
SDPS induces subgap resonances in the charge conduc- 
tance of the quantum point contact. For high transparen- 
cies, these resonances are smoothed, but the shape of the 
signals remains extremely sensitive to SDPS. In Ref. [s^l 
it has been found that when F is strongly polarized, the 
peak in the conductance of the F/I/s- wave superconduc- 
tor contact can be restored at a zero potential. 

Such strong influence of SPDS on transport properties 
of hybrid structures with ferromagnetic elements allows 
one to suppose that they may be successfully used in 
experiments on Andreev spectroscopy of ferromagnets, 
superconductors and in various applications in the field 
of nanospintronics. 

This paper is devoted to a theoretical study of the 
SDPS influence on the Andreev reflection and charge con- 
ductance of a point F/I/d-wave superconductor contact. 

Superconductors with the d-wave symmetry (the 
dx^-y^ symmetry of the order parameter is considered) 
have an internal, momentum-dependent phase, which 



2 



strongly influences the transport properties of contacts 
between them and other materials. In Ref. i34] it was 
shown that when the angle 7 between the a axis of a su- 
perconducting crystal and the normal to the surface of 
the high-ohm interface is 7r/4 (the {110} orientation of 
the d-wave superconductor) , then a bound state is formed 
on the Fermi level near the high-ohm interface. This zero- 
energy bound state resulting from the repeated Andreev 
reflections [H, [s^ causes a sharp peak at a zero poten- 
tial in the dependence of the charge conductance of the 
N/I/d-wave superconductor on the applied potential [st! ]. 

The first theoretical study of spin-polarized tunnel- 
ing spectroscopy of F/I/d-wave superconductor junctions 
was performed in Refs. [ssl. [sol [40j. It has been found 
that the subgap charge conductance behavior is quali- 
tatively different from a nonmagnetic case. In particu- 
lar, it has been found that for the {110} orientation of 
the d-wave superconductor the zero-potential peak in the 
charge conductance is suppressed by the exchange inter- 
action due to the suppression of Andreev reflections and 
that it splits into two peaks under the influence of the 
exchange interaction in the insulator. 

The influence of SDPS 0f and 9^^ associated to the 
electronic reflection and transmission amplitudes r and d 
on the contact on the charge conductance of the F/I/d- 
wave superconductor contact in Refs. [SSj, |39j, |40j is not 
studied. 

The main result of this paper is that SDPS can sup- 
press the zero-potential peak and restore finite-potential 
peaks in the charge conductance of the F/I/d-wave super- 
conductor contact for the {110} orientation of the d-wave 
superconductor and, on the contrary, can restore the 
zero-potential peak and suppress finite-potential peaks 
for the {100} orientation of the d-wave superconductor. 
This takes place because due to the interference of one 
part of trajectories of electron-like and hole-like quasi- 
particles reflected by the pair potential and the interface, 
spin-dependent bound states are formed near the Fermi 
level, whereas due to the interference of the other part of 
trajectories spin-dependent bound states are formed in 
the vicinity of edges of the energy gap. Spin-dependent 
amplitudes of the Andreev reflection probability and en- 
ergy levels of spin-dependent Andreev bound states are 
also found. 

This work illustrates that the study of the influence 
of SDPS on the charge conductance of the point F/I/d- 
wave superconductor contact can provide an interesting 
insight in the spin-dependent transport. 

A theoretical possibility to study the influence of SDPS 
associated to the electronic reflection and transmission 
amplitudes acquired by electrons upon scattering on the 
potential barrier on the I —V characteristics of supercon- 
ducting weak links with ferromagnetic elements appeared 
after the boundary conditions (BCs) for the quasiclassi- 
cal GF were obtained. In Ref. [4l|, BCs for the quasi- 
classical GF for two metals in contact via a magnetically 
active interface in terms of an interface scattering ma- 
trix were derived. In Ref. |29|, BCs for the retarded and 



advanced quasiclassical GFs were obtained in terms of 
Riccati amplitudes [43,|43|]. In Ref. [31,], BCs in terms of 
Riccati amplitudes were obtained for the nonequilibrium 
quasiclassical GF. In Ref. 'ii'l, quasiclassical equations of 
superconductivity for metals with a spin-split conduction 
band were derived and BCs for the temperature quasi- 
classical GF for the F/S interface were obtained. The 
model interface was the same as in Refs. |45,] . 

In this paper, calculations are carried out using quasi- 
classical GFs and the relevant BCs obtained in Ref. Hi- 



ll. FINDING DIFFERENTIAL CONDUCTANCE 
OF A POINT FIS CONTACT 

A. The general expression for differential 
conductance of a point contact through 
quasiclassical GF 

In hybrid F/S structures the Andreev reflection is mod- 
ified. The reflected hole has some parameters (for exam- 
ple, the velocity modulus and the phase shift) different 
from those of the incident electron because it moves in a 
subband with an opposite spin. Such spin-discriminating 
processes due to the exchange fleld in a ferromagnet lead 
to the formation of spin-dependent Andreev bound states 
inside the gap [1^, |33|. As a result, the spectral density 
of the charge conductance Gpis of the FIS contact at 
a zero potential is no longer a symmetrical function of 
energy e. The condition of the time reversal invariance 
has the form Gpisis, a) — Gfis{— e, — a). The general- 
ization of the charge conductance expression [i^ for this 
case results in the following formula for Gfis{V) [33j : 



GpisiV) = 



327r2T 



dPlI 



(2 Try 



de X 



coth^( 



s 'zy s • z 



2T I 



(1) 



In Eq. (Ill), F is the potential; A is the contact area; e 
is the electron charge; T is the temperature; fz is the 
Pauli matrix; p|| is the momentum in the contact plane; 

{(js-, Ts) and {cja, Ta) are quasiclassical retarded (R) 
and advanced (A) GFs symmetric (s) and antisymmetric 
(a) [33*1 with respect to the projection of the momentum 
p on the Fermi surface on the x axis, being perpendicu- 
lar to the contact plane, composed according to the rule 
f,(,) = 1/2 [T(p,) ± f{-p,)]. 

Besides the matrix quasiclassical GF g, which equa- 
tion is analogous to that derived in Ref. [11], equation 
(H]) includes the matrix GF T, describing the interference 
of waves incoming to the interface and outgoing from it. 
The function relation with the matrix one-particle tem- 
perature GF and equations, which the function obeys. 
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are presented in appendix. Calculations in Eq.((T]) are to 
be carried out on the boundary of any contacting metal. 



B. Finding quasiclassical GF 

Let us assume that the barrier with the width d is 
located in the region —d/2 < x < d/2, the supercon- 
ductor occupies the region x > d/2, and the ferromagnet 
occupies the region x < — d/2. To find GFs for each 
metal, one has to solve quasiclassical equations of su- 
perconductivity for metals with a spin-split conducti vity 
band simultaneously with their BCs derived in Ref. [33] : 



ox 

sign(p^)^t -I- 



1 5 i„ 



5«x 



[K, g]- = 0, 



2 ' 'dp' 
_i - - _i 

K ^ ~iVx^ {iSnTz + A - ^)Vx'' - iiPx - TxPxTx)/'^, 

A = A{x, p), [a, b]± ^ ab ± ba. (2) 

In this section, £„ — {2n + 1)ttT is the Matsubara fre- 
quency; fx and fz are the Pauli matrices; p = {x, y) 
are coordinates in the contact plane; S is the self-energy 
part; g are matrix temperature GFs: 
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Moreover 
A = 



9 a fa, — Q 
f—a,a 9— a,— a 



, 9 



9> 
9< 



A(x, p) 
-A*{x, p) 



Px 



Px,c 





Px > 0, 
Px <0 





Px,-c 



where A{x, p) is the order parameter; Px^a and p|| are 
projections of the momentum on the Fermi surface on the 
X axis and the contact plane, respectively; Vx =Px/m 
and = p\\/m. 

BCs for a specular reflection of electrons from the 
boundary p|| — sin-i^i — p^ sinz^i = ps sini?5 have the 
form 33]: 



(If )<i = i9a)d, (Tf = (T-),, 



(3) 



= \[k'o,{gt)d + ai{gt)n, 
-^t = {Rc.yHk)d+a2{97)n, 



where = 1/2 [.g„(^) ± ]. Functions t^^^^ are 

determined analogously. In Eq. ([3]) and below, the in- 
dices d and n denote the diagonal and the nondiagonal 
parts of the matrix {T^(n) = 1/2 JT± r^Tr^]), respec- 
tively. GFs g are connected with GFs, being the 



solutions of Eq. ([2]), by the following relationships [33j : 



(5f)« = i9f)n cos{Oa.)+ifz{g!)n sin(0„) 
i9a)n = (fff)n cos( „ ) + if, (.gf )„ sin(0„) 



(.9f)n = (5f)n cos(/3;) +if,(5f)„ sin(/3; 
(5f)n = i9a)n cos(/3;) -f if, (gf )„ sin(/?;;) 



r /I r 

d nd \. or Q '^ — a 
a~^~a)i Pa- 7, 
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The diagonal parts of matrices g are equal to the corre- 
sponding matrices g. An explicit form of other functions 

T is not needed. These functions are found from BCs. 
Coefficients a, are: 



"1(2) = 



1 + 



«3(4) = 1 - ^R^Ri ± V^T-Di )■ 

When solving Eqs. jH, let us assume that the order pa- 
rameter is a step function, being zero in the ferromagnet 
and finite in the superconductor. Then for S metal the 
solution is as follows: 

g(a;,p) = e-"'sn(j5.)K(x-|)(j,(p)gSig„(p.);$-(x-|) _^ ^^^^p^ 

T(a;,p) = e-'5ign(j5^)K(a;-|)Yg-sign(35^)it(a;-^) 

t = t(x = 0,p). (5) 

Matrices C'o(p) are the values of GFs g far from the F/S 
boundary: 



A(P) 



/ 



tA*ip) 



^*A(P) 

-£n 



5/ 

Ad(r) cos(2i?s - 



-\MpW 

27). 



(6) 



In Eq. ([6]) Arf(r) is the maximum value of the order 
parameter at temperature T; i!) s is the angle between 
the electron momentum in the superconductor and the 
X axis, being perpendicular to the contact plane, and 7 
is the angle between the crystal a axis of the d-wave 
superconductor and the x axis. 

For F metal the solution has the same form as Eq. ^ 
except for changing the exponent argument from (a; — 
d/2) to [x + d/2); = fz £„/]£„]. 

GFs 5'^ in Eq. ^ have to tend to at a; ^ -|-oo 
and GFs g^ to Cff (p) at a; ^ — 00. By matrix multipli- 
cation in Eq. ([5]) and in corresponding equation for 5^, 
we find that for the above to hold it is necessary that at 
X — -\- d/2 and at x — — d/2 the relationships 

Co^(p)C^(p) =(7^(p)(7o^(p)=sign(p,)C^(p) (7) 
Co^(p)C'^(p) =~(7^(p)(7o^(p)= -sign(p,)C^(p) 
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are fulfilled respectively. It follows from these relation- 
ships that 

9a = C'f +(^0,a^ 9a =C'f, (8) 

where 

^ - (1 + Cla)iCls)'\ X = f, {X)a + (1)„. 



In Eq. 



C 



0, s(a) 



are symmetric and antisymmet- 



ric combinations of the matrix C'o^(p) with respect to 
the projection of the Fermi momentum on the x axis: 
Cls(a) = ^mClip.) ± Cl{-p,)], X ^l. Matri- 
ces T"^^^) satisfy the relationships: 

Cff (p)t^(p) = t^(p)C'o^(p) = -sign(p,)(p)t^(p) 
(7o^(p)t^(p) = t«(p)Co^(p) = sign(p,)(p)t^(p), (9) 

being the condition for the hmctions T''^(a;,p) and 
"T^ [x, p) to tend to zero when x tends to — oo and -I- oo, 
respectively. It follows from Eq. ^ that the func- 
tion (T(x)^)„ ~ const — 0, because for a ferromagnet 
[K, (T)„]+ = and 'T^{x, p) has to tend to zero when x 
tends to — oo. Then from the BCs ^ and relationships 
(gl) it follows that: 



0^3(9 a )n = a4(.ga )r 



= "2(5^ )r 



(10) 



From the first equality in Eq. (jlOp we find the relation 
between functions (g^)n and (gf )«: 



(ID 



^ (If) 



By substituting this relation into the second equality in 
Eq. PH)) and using the relations dU and ([5]) we find 



9a 9a ~ ^ 



Z = {1- ^R^i) [Xd cos{e^) + 1 sin((?„)] (11) 
+ (1 + ^/R^Ri) sign(e„) [cos(6l„) + i Xd sin(6'„)]. 

Knowing (g^)n, from Eqs. ([3]) and ^ we find functions 
gf,Tf and T ^ necessary for calculation of the conduc- 
tivity in the Eq. (P) and calculate the conductance at 
the ferromagnet side. 



C. Differential conductance of a point FIS contact 

After carrying out the analytical continuation in func- 
tions {g^)n,gs ,Tf (substitution i e„ for e±S for 
retarded and advanced GFs, respectively), we obtain an 



expression for the charge conductance a pfs{V). For an- 
gles 7 = and 7 = 7r/4 crp^s{V) is as follows: 




cosh^(^^) 

|A(<»s)| 

de 

2r 



cosh2(£±|^ 



+ 
(12) 



DfD 



1^1 



cosh^(^) cosh'(^^) 



|A(^5)|^ 



For 7 = 0: 

A(7?5) = |A<i|cos(2i95) (13) 
= [£(1 -W)+ e(l + W)r+ AW |A(i?s)|' sin^e^) 
Zi^ = [l + 2M/^cos(26»„) + W^]\A{iSs)\'^ - 4We^cos{2ea) 

16W^{\A{ds)\^ -£^)s^ sin^(2 0^) 

[1 + 2W cos(2 6ia) + W'^] I A(i?s)P - ^We^ cos(20„) 
W = /R^?!; ^ = - \A{dsW- 

For 7 = 7r/4: 

A(79s) - |A<i|sin(279s) (14) 
Z^y = [e(l + VK) + e(l -W)^- AW \A{ds)? sin2(0„) 
Zjl = [1 - 2Wcos{26a) + W'^]\A{ds)\'^ +4We^cos{2ea) 

161^^ (I A(t^g)P - £2)£^ sin'(2 0^) 

[1 - 21^ cos(2 6ic«) + 1^2] I A(z95)|2 -|- 4M^£2 cos(26'„) ' 

For 7 = 0, when 9^ — 0, the expression for the con- 
ductance obtained in Ref. [s^l follows from Eq. p2|l . 
In the case of a nonmagnetic metal, when D-^ = this 
expression is the same as that obtained in Ref. [i^ , and 
for D = 1/(1-1- Z"^) this expression is the same as that 
obtained in Ref. For 7 = 7r/4, when 6a = 0, 

the expression for the conductance obtained in Ref. [47j 
follows from Eq. p2)) . 



III. ANDREEV REFLECTION 

The calculation of quasiclassical GFs in the expression 
for the conductance allows one to conclude that for en- 
ergies lower than |A(t95)| {e'^ < |A|2), the following 
relation is true: 



^T^T,T^:T,]=A[-g^Trg^T^]^l. (15) 

The comparison of the form of under-gap conductances 
in Eq.(IT]) and that of the corresponding Eq. (25) in Ref. 
[46i] shows that the matrix elements (.g^)^ and (5^)^ 
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are the amplitudes of the Andreev reflection probability 
a{e, 9a) in FIS contacts. Let us take the matrix elements 
of (ff^)^ given by Eq. ^ as a{e, d a): 



a(7, e, 9a) = 



Z{l) 



(16) 



where 



Z(0) 



(1- yiV^)[£COs( 



f z(i+v/i!p?r)[vw^ 

Z(^/4) = 



V|A(i?s)P-£2sin(0„)] 
■2 cos(0Q,) +e si^(^^Q,)]. 



(1 + /R;Rr)[£Cos(0„) - v/|A('?s)p-£'sin(0„)] 
+ i (1 - /Rp?r)[V|A(^9s)P - £2 cos(0„) +e sin(0c«)]. 

The presence of the imaginary part in functions 
0(7, £, 9a) means that Andreev reflection is accompa- 
nied by the phase shift. The Andreev reflection proba- 
bility ^^(7, e) {Aa{^,e) = a(7, s, 6a)a*{^, e, 9a)) is: 



Aa{l,e) 



DlWs)[' 
\Z{l)? 



(17) 



|Z(0)p = [l-v^]'|A(^^,)|^ 



+ A^R^^[./\K{ds)\ 
|Z(7r/4)p 



■2 COs{9 a) 

1 2 



e 

[1-/rP?J^|A(795 



e sin( ( 



+ A^fR^^W\^{ds)\^ 



sin(6'Q) — e cos(0q) 



From this equation it follows that: (1) the spin- mixing 
angle 8 used in Refs. [28', '29*1 corresponds, in our nota- 
tions, to 6i„ (for S/F/S and N/F/S contacts Q = e'^-9l^ 

9'^ -9f di, [3^); (2) for 7 = 0, when 9a < the 
Andreev reflection probability of the electron excitation 
with the spin projection a is larger than that of the hole 
excitation; when 9a > 0, the Andreev reflection proba- 
bility of the hole excitation with the spin projection a is 
larger than that of the electron excitation; for 7 = 7r/4, 
the relation is reverse; (3) the Andreev reflection proba- 
bility has maxima at e = £^ (7) (at values of the energy 
of electron (hole) excitations corresponding to the energy 
levels of spin-dependent Andreev surface bound states). 
The energy of spin-dependent bound states is: 



e'^(-) 



- \A{ds)\cos{9a) for {tt/2) >9a>0, 



\A{^s)\cosi9a) for - (V2) < < 0, 
|A(i?s)||sin(0a)| for (7r/2) > 0„ > 
|A(7?s)||sin(0„)| for - (7r/2) < 0„ < 0. 



(18) 



Spin-dependent Andreev surface bound states are 
formed in a superconductor due to the interference of 
electron-like and hole-like particles with different SDPS. 
One may demonstrate it by using a phenomenological ar- 
gument in Ref. [36j . Let us consider diagrams in Fig. 1 
corresponding to Andreev reflection of an electron with 




FIG. 1: Structure of the diagrams corresponding to Andreev 
reflection in the superconductor: diagram a), one-act pro- 
cess; diagram b), two-act process. The vertex Q is Andreev 
reflection of electron-like (solid lines) and hole-like (broken 
lines) quasiparticles by the pair patential. The vertex • is the 
normal reflection of electron-like and hole-like quasiparticles 
by the barrier potential. When the solid line transforms into 
the broken line, Q denotes the vertex Pa^-a- When the bro- 
ken line transforms into the solid line, O denotes the vertex 
P-a,a- Parameters da, da, and fa are related as follows: 

da ~ d a P X / P X a] ^ ~ 1^ a ^ / ^ a] Da ~ dada [45| ] . 



the spin projection a and the energy less than |A| trans- 
mitted from a ferromagnet into a superconductor. The 
analysis of these diagrams and their summation makes 
it possible to obtain the following expression for a phe- 
nomenological expression of the amplitudes of the An- 
dreev reflection probability a(e, 9a)'- 

a{e, 9a) = darfl„/3„^'ijl + rl„ f „ 
+ (r_f„/?„^^„/?„^%)^ + ...] = 



'—a'oih'a,—at^—a,a 
DaD-a Px a/P X - a ^ ^ " Pa. -a 

(19) 



e''«= - e-^^-^/Rji~al3-^_aPh^^a 
The corresponding probability of Andreev reflection is: 

/) ^ _ ^ci^~a P x,a/ P X, - a ya.-a Ha. -a 
Q = ^/R^a [C0S(2 9a)[(i^!'.a P-a. a + /^^f-a P-ata] 

+ ^ sin(2 9a)[P*a:'\ P-^U " P-l al] (20) 

By comparing formulas (|161ll7p . derived using quasiclas- 
sical GFs, with formulas obtained using phe- 

nomenological arguments, we find the expressions for the 
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vertices and So for 7 = 7r/4: 




|A(^?5)| 

^ A*(z9s) 

|A(t?5)| 



(21) 



For 7 = the expression for the vertex /Jj^^J ^ is of an 
opposite sign. It foUows from formulas (|^D|) and (pij) that 
in the absence of the interferential term Q the probabil- 
ity of Andreev reflection is a constant (independent of 
the energy e) quantity. The interference of electron-like 
and hole-like particles reflected by the pair potential and 
the interface results in the formation of spin-dependent 
Andreev surface bound states. For 7 = at = the 
maximum in the probability of Andreev reflection is at 
s = ±|Ad| as in [46]. At 9a = ±7r/2 spin-dependent 
Andreev surface bound states with the width F: 



(l-^/i?^RI)|A(^s)| 



(22) 



0.5 

0.0 h 
-0.5 
-1.0 
-1.5 



k/p,= 0.85 



p,d = 2 



5=0.85 
0.65 




0.0 0.2 0.4 0.6 0.8 1.0 
C0S(i3^) 



are formed at e = on the Fermi level. For 7 ~ 7r/4 the 
spin degeneracy of the level on the Fermi surface |34!| at 
9a is removed. Two energy levels symmetric with 
respect to the Fermi level are formed inside the energy 
gap. 



IV. APPEARANCE OF ANDREEV BOUND 
STATES IN CONDUCTANCE OF THE FIS 
CONTACT 



We present below the results of numerical calculations 
of the charge conductance of the FIS contact taking into 
account the phase shifts. In the numerical calculations 
the relation between Fermi momenta of contacting met- 
als was the following: ps — (pt + pi)/2. Calculations 
are carried out for a rectangular barrier with a height U 
counted from the bottom of the conduction band of a su- 
perconductor. The electron wave function in the isolator 
x(,x) is as follows: 



Xix) = Ci exp(^a;x) -I- C2 exp{-n^x), 



where ^ix ^ yk'^ + pf, k'^ = '2mi,{U - E§,) ; E§, is the 
Fermi energy of a superconductor, rrib is the mass of an 
electron in a barrier. In this case the expressions for 9'^ 



FIG. 2: Dependence of the angle on cos(i?|) for various 
values of the polarization of a ferromagnet 5 ((5 = Pi/pi < 1). 



and 9 a have the following form: 



^ d 



ad , Apf.aP^ -Mx)tanh(/i:,d)\ 

9a = arctan ^ — —- (23) 

V {pi^a + P%) j 

-Xr , /2^:,pf,„ ^ (pf)2]tanh(/x^d)\ 
9^ ^ arctan I ■ I 

-^[^it - (pf pi^„)2]tanh2(Ai,d), 

so that the angle 9 a^ [9a = (6'a-6'-a)/2-(6'a-6'^ a)] = 
{Oa-SLa)/'^-i^i-^-a) ^ocs uot depend on the location 
of the barrier. 

Figure 2 shows dependences of the angle 9^ on 
cos(?9|). All angles are connected by a specular reflec- 
tion p|| = p|sint9| = j3|sin?9| = pssin'ds- Figure 2 
shows that the angle 9 1 , being a combination of phase 
shifts 9^ and 6*^, has a jump for a part of electron tra- 
jectories transmitted through the contact region. The 
jump in the dependence of the angle 6'^ on cos(i?x) is 
due to the jump in the dependence of the phase shift 9^ 
on cos( -d i ). 

By setting p^^^ ^ Px,i ^ Px = Px Eqs.(1221), we 
get the following expression for phase shifts 9^ and 9"^, 



7 



which correspond to the N/I/S contact with the same 
t:?„„™; — 4-„ „N _ „s normal metal {p^) and 

:p^cos(z9)): 



Fermi momenta p 
superconductor {p^ 



p' 



arctan 



N 



9 = arctan 



[/x2-(p^)2]tanh(/x,d) 
(pf)2]tanh(Ai,d) ~, ; 



(24) 



2nxP 



N 



One may see from this equation that the phase shift 9'^ 
is continuous for the trajectory with p^ = Hx: and the 
phase shift 0'' differs by tt for electron trajectories with 
p ^ , being somewhat larg er /i X and somewhat smaller /i ^ • 
For the N/I/S contact, phase shifts 9^ and have jumps 

for the same trajectory /u a;- As a result, the angle = 0. 
For the F/I/S contact, the trajectories, at which phase 
shifts 9^ and 9^^ have jumps, are different. These critical 
trajectories are the solutions of the transcendent equation 
Za = (see Eq. (23)). An analysis of the numerical 
solution of this equation allows us to state that if fed < 2 
and k/pi < 1 there is always a set of trajectories, for 
which the phase shift 9^ has a jump and the phase shift 

9^ has no jump, or vice versa. Thus, the angle 9a has 
a jump being equal to ±7r/2. Figure 2 shows the case 
when the phase shift 9^ has a jump and the phase shift 
0| has no jump. 

As the polarization of a ferromagnet 6 (6 — Pi/p] < 
1) increases, the part of electron trajectories with phase 
shifts experiencing a jump increases as well, and at high 
values of the polarization of the ferromagnet for all elec- 
tron trajectories 6*1 ~ — 7r/2,0| ~ +Tr/2. For a rectan- 
gular model of the potential barrier, the angle 9 a is of 
an order of (=F7r/2) only for kd < 2 and k/p^ < 1. At 
fc/pf > 1, the angle 9a << 1- 

Figure 3 shows the results of numerical calcula- 
tions of the normalized conductance of the FIS contact 
o'F/s{V)/ao for the {100} - oriented d - wave supercon- 
ductor taking into account and not taking into account 
the phase shifts. Not taking into account the angle 9 a 
(dashed lines) , the plots illustrate the suppression of An- 
dreev reflection due to a decrease of the number of con- 
ducting channels determined by the number of conduct- 
ing channels in the subband with a lower value of the 
Fermi momentum (in this case it is p |) as the polarization 
of the ferromagnet increases. The appearance of electron 
trajectories with a jump of the angle 9 a forms Andreev 
bound states on the Fermi surface (£^(0) = 0) (fT5|) . It 
results in the appearance of a peak at a zero potential in 
the dependence a p^si^)- As 6 increases, the part of elec- 
tron trajectories around the normal to the contact plane 
participating in the formation of levels close to the Fermi 
level of a superconductor increases (Fig. 1). However, 
it does not compensate the decrease of the conductance 
at a zero potential due to the decrease of the number of 
conducting channels. As a result, with increasing polar- 
ization of the ferromagnet the zero-potential peak in the 
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FIG. 3: Dependence of the normalized conductance 
ap/s{V)/ao (cro = e^P|/7r^) from Eqs.([l2l) and ([13} as a 
function of the applied potential for the {100}-oriented d - 
wave superconductor (7 — 0) for various values of the polar- 
ization of a ferromagnet 5 not taking into account (dashed 
lines) and taking into account (solid lines) the phase shifts. 



dependence cr f/s{V)/'^o decreases. 

Figure 4 shows the results of numerical calcula- 
tions of the normalized conductance of the FIS contact 
a p/s{V)/o'o for the {110} oriented d - wave superconduc- 
tor (7 = 7r/4). A part of electron trajectories without a 
jump of the angle 9 a forms the Andreev bound state on 
the Fermi surface manifested in the zero-potential con- 
ductance peak. The other part of electron trajectories 
with a jump of the angle 9 a forms the Andreev surface 
bound state with an energy of about jA^ sin(2i?5)| man- 
ifested in the conductance peak at the potential close to 
\Ad\. At increasing polarization of the ferromagnet all 
electron trajectories have a jump of the angle 9a- As a 
result, the conductance peak at the zero potential disap- 
pears and that at the potential close to ±\Ad sin(2'!95)| 
increases. Plots in Fig. 4 demonstrate a tendency of the 
conductance peak to the decrease at the zero potential 
(the decrease of the part of electron trajectories forming 
the level on the Fermi surface) and the increase of the 
conductance at the potential close to the edge of the su- 
perconducting gap (the increase of the part of electron 
trajectories forming the jump of the phase shift 9a) with 
increasing polarization of the ferromagnet. 
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FIG. 4: Dependence of the normalized conductance 
<^ F/siy)/'^o from Eqs.((T2} and (|14|) as a function of the ap- 
phed potential for the {110}-oriented d - wave superconductor 
(7 = 7r/4) for various values of the polarization of a ferromag- 
net 5 not taking into account (dashed lines) and taking into 
account (solid lines) the angle 9-\. 



V. CONCLUSION 

In this paper, the influence of SDPS associated to 
the electronic reflection and transmission amplitudes ac- 
quired by electrons upon scattering on the potential bar- 
rier on the Andreev reflection probability of electron and 
hole excitations for a ferromagnet/isolator/d- wave super- 
conductor contact and the charge conductance of the 
ferromagnet/isolator/d- wave superconductor contact as 
a function of the applied potential have been studied. 
Spin-dependent Andreev bound states in a superconduc- 
tor are found. It is found that for parameters of a po- 
tential barrier kd < 2 and k/pij < 1 there are always 
two groups of trajectories of electron-like and hole-like 
quasiparticles, such that due to the interference of one 
group of trajectories of electron-like and hole-like quasi- 
particles reflected by the pair potential and the interface, 
spin-dependent bound states are formed near the Fermi 
level, whereas due to the interference of the other group 
of trajectories spin-dependent bound states are formed in 
the vicinity of the edges of the energy gap. As a result, 
SDPS can suppress the zero-potential peak and restore 
finite-potential peaks in the charge conductance of the 
F/I/d-wave superconductor contact for the {110} orien- 



tation of the d-wave superconductor and, on the contrary, 
can restore the zero-potential peak and suppress finite- 
potential peaks for the {100} orientation of the d-wave 
superconductor. The fitting of Eq. to the experi- 

mental dependence of the charge conductance of the FIS 
contact on the applied potential makes it possible to de- 
termine the polarization of a ferromagnet. 
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APPENDIX: DETERMINING QUASICLASSIC 
GFs g AND t. DERIVING EQUATION (2). 

Let us start with equations for equilibrium thermody- 
namic GFs in the matrix form [48j, taking into account 
the spin splitting of the conduction band: 

1 - A . 

lenTz + — — -I- A(r) + A - S G(e„, r, r') 



2m dr^ 
S{r-r'). 



(A.l) 



Here S is the self-energy part which includes the scatter- 
ing by non-magnetic impurities and phonons ^48| . An 
explicit form of this term is not needed for deriving the 
quasiclassic equations. G(e„,r, r') is the matrix temper- 
ature GF: 

,2 



F 

G — (y — r 



1 
2m 



Pa 



P- 



fz is the Pauli matrix; e„ = {2n+ l)7rT is the Matsubara 
frequency, a is the spin index; A(r) is the order parame- 
ter (as defined below equation (2)); pa is the Fermi mo- 
mentum; m is the electron mass; r = (x,R), R = (y, z); 
X-axis is perpendicular to the contact plane. 

Passing to coordinates p and p {p — p—p', 2p = R+R') 
in Eq. (|A.1|) and performing Fourier representation with 
respect to the p coordinate, the following equation for 
G{x,x') = G{x,x' , p,p\\,en) {p\\ is the momentum in the 
contact plane) is obtained: 



1 92 



d 



2m dx^ 2 dp 
— S{x — x'), 



Px 

2m 



(A.2) 



21 1/2 



In Eq. (A2): V|| = P||/m, p^ = [pi - p\\' 

Then the Zaitsev representation generalized for the de- 
scription of metals with a spin-split conduction band is 
used for the function G{x^x'): 

2 

G(x,x')= Mx)Gn,m{x,x')Al{x'), (A.3) 



n, m— 1 



Px 



_ (Px,a 

Px.-a, 
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Representation (|A.3|) explicitly takes account for oscil- 
lating terms present in the function G{x,x') and waves 
of the exp[±i (p|X + pj.x')] type, arising from partial re- 
flection of the first electron of the superconducting pair 
from the interface |49|]. Functions G n,m{x,x') change at 
distances of an order of the mean free path of electrons 
in a metal. By substituting Eq. (jA.3p to Eq. (|A.2|1 and, 
neglecting the second a;-derivative, we obtain an equation 
for slow changing functions G k.n{x,x'): 

Ak{x) (iSnT, - ii-l)''v.^-^ + + Mx) (A.4) 



Px 



' dx 2 dp 

-t) Gkn{x,x')An{x')=5{x-x') 

/ m 
Analogously, an equation conjugate to (|A.1[) gives: 

Ak{x)Gkn{x,x')An{x') ^i£„T^ + + H^' ) - S 

+ Ak{x)^^^^^^^An{x')i{-irv^ = 5{x - x'). (A.5) 

In Eqs. (A4) and (A5) let us pass to functions go = 
go{x,x') = gQ{x,x' ,p,p\\,en) and Tq = to(x,x') = 
To(a;, a;', being continuous at a point x = a;', 

by using formulas: 



ffo = 



,95" = 2i^/vxGii{x,x')^/vx - sign(a; - x'); p^, > 



g^ = 2iyJv^G22{x,x')^/vx + sign[x - x'); Px < 

2iy^Gi2{x,x')y^, Px > 

(A. 6) 



T> 



[ T;5^ = 2iy^VxG2l{x,x')y^Vx, < 0. 

Let us call the obtained equations as (A4') and (A5'), 
respectively. By subtracting equations (A5') from equa- 



tions (A4') when n = k and adding equations (A4') and 
(A5') when n ^ k, one may get equations for functions 
go{x, x') and To(a;, x'). In these equations we set x — x' . 
Finally, the following equations are obtained: 



sign(p.)S(x)^ B*{x) + ^ ^^Bix)[v-\ go(x)] + B*(x) 



vii d 



-[ko, B{x)goB*{x)]- = 0, 



(A.7) 



sign{px)B{x)^B{x) + ^^^B{x)[v-\ to]- B{x) 
+ [ko,B{x)foB{x)]+=0, (A.8) 

B{x)^e'''s^^P^'>P^^, ko^ -z{>;'(i£„f, + A-E){};^ 
[a, b]± ^ ab±ba. (A.9) 

Considering that in expression for B{x) the matrix 
can be written with the help of the Pauli matrix Tx as a 
sum of two components proportional to the unit matrix 
and Pauli matrix f^: 

Px = {Px + TxPxTx)/2 + {Px - TxPxTx)!'^, (A. 10) 

and going in Eqs. (A7) and (A8) to functions g {g = 
g{<^n,P\\,P,x)) and t (t EE t(e„,p||,p,x)) by formulas: 



g = g**5'Sn(Pm)f2 2;^^g-isign(p^)r2a: 
Y = gisign(pi)na:Yggisign(p^)f2a; 
^ P x ^xP x^x 



(A.ll) 



one obtains equations 

If quasiclassic GFs g and T are independent of the p 
coordinate, the condition g^ — 1 is met. 
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